Damped Bogoliubov excitations of a condensate interacting with a static thermal 

cloud 
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We calculate the damping of condensate collective excitations at finite temperatures arising from 
the lack of equilibrium between the condensate and thermal atoms. Since a static thermal cloud 
already produces damping in our model, we ignore the non-condensate dynamics. We derive a set of 
generalized Bogoliubov equations for finite temperatures that contain an explicit damping term due 
to coUisional exchange of atoms between the two components. We have numerically solved these 
Bogoliubov equations to obtain the temperature dependence of the damping of the condensate modes 
in a harmonic trap. We compare these results with our recent work based on the Thomas-Fermi 
approximation . 

PACS numbers: 03.75.Fi. 05.30.Jp 67.40.Db 



I. INTRODUCTION 

In this paper, we calculate the damping at finite tem- 
peratures of collisionless condensate collective modes due 
to interactions with a static thermal cloud. Our treat- 
ment starts with the general theory developed in Ref. ||] 
and is an extension of the static Hartree-Fock-Bogoliubov 
Popov (HFB-Popov) approximation |g, ^, Q] to include 
the effect of the lack of diffusive equilibrium between 
the condensate and non-condensate. The static thermal 
cloud acts as a reservoir that can exchange particles with 
the condensate so as to bring the condensate back into 
equilibrium, resulting in a damping of the collective os- 
cillations of the condensate. We do not expect, in gen- 
eral, that the coUisional damping we find from a static 
thermal cloud will be much modified when we treat the 
non-condensate dynamically. Because we are treating the 
non-condensate statically. Landau and Beliaev damping 
are not contained within in our description. 

Our present work is the natural extension of a re- 
cent study H of the damping of collisionless conden- 
sate collective modes at finite temperatures. This was 
based on a simplified model treating the time-dependent 
Gross-Pitaevskii equation for the condensate within the 
Thomas- Fermi approximation (TFA). We described the 
damping in terms of the density 6nc(r,t) and velocity 
Svc{v,t) fluctuations of the condensate oscillations. At 
T — 0, where the non-condensate component is absent, 
Stringari [pj showed that a simple analytic solution of the 
condensate collective modes could be obtained within the 
Thomas- Fermi approximation. In Ref. [pi, we formally 
generalized these "Stringari" solutions to finite tempera- 
tures by assuming that the mean field of the static ther- 
mal cloud has a negligible effect on the condensate oscilla- 
tions. We then calculated the damping due to the lack of 
equilibrium between the condensate and non-condensate 
to first order in perturbation theory, using the Stringari 
modes as a zeroth-ordcr solution. In Ref. B we compared 
this damping directly to Landau damping and found it 
gives a 30% to 50% correction to the total damping for 
typical experimental parameters. 



The main focus of the present paper is to generalize 
our earlier work ||5|] to include the same inter-component 
coUisional damping mechanism directly in the Bogoli- 
ubov equations, still treating the non-condensate stati- 
cally. This allows us to check the range of validity of the 
TFA made in Ref. |). The TFA is valid when the num- 
ber of condensate atoms Nc is large, that is, when the 
condensate radius Rtf is much larger than the harmonic 
oscillator length tho 



tw/Rtf = (15iV,a/rHo)"'/' « 1, 



(1) 



where a is the s-wave scattering length. The TFA will 
break down as one approaches the BEG transition tem- 
perature Tbec where Nc becomes small. It is in this 
temperature region where we must solve the coupled Bo- 
goliubov equations for our model, as we do in this paper. 
From the time-dependent GP equation describing the 
fluctuations of the condensate order parameter, i54>(r,i), 
we derive generalized Bogoliubov equations for the col- 
lective mode amplitudes Uj (r) and Vj (r) and complex en- 
ergies Ej , which we solve numerically and compare to the 
predictions of our earlier TFA calculation [p|. We also 
give a careful discussion of the formal properties of cou- 
pled Bogoliubov equations which include damping, since 
this may be of more general interest than our specific 
model calculation. 



II. STATIC POPOV APPROXIMATION 

Our starting point is the finite T generalized GP equa- 
tion derived in Ref. pi (see also Refs. Iq, H].) 



ih 



a$(r,i) 
dt 



2rn 
iR[v,t) $(r,i) 



V^ + C/cxt(r) + .g?^c(r, t) + 2.gn(r, t) 



(2) 



where the interaction parameter g — AttH a/m, a is the s- 
wave scattering length, nc(r, t) — |$(r, i)p, and h{r, t) — 
J dp/(p, r, t)/{2Trh,)^ is the non-condensate local density. 



The term i?(r, t) in (g) describes the exchange of atoms 
between the condensate and normal gas and is given by 



Rir,t) 



2nc{v,t) J {2TTh) 



dp 



:Cu[f{p,r,t),^r,t)]. (3) 



This involves the collision integral C12 [/, $] describ- 
ing collisions of condensate atoms with the thermal 
atoms, which also enters the approximate semi-classical 
kinetic equation for the single-particle distribution func- 
tion (valid for fcer > Mco) 

^%^ + ^•Vr/(p,r,t)-Vt/(r).Vp/(p,r,i) 
mm 



= Ci2[/,$]+C22[/]. 



(4) 



Here the collision integral denoted by C22[/] describes 
binary collisions between non-condensate atoms. It does 
not change the number of condensate atoms and hence 
does not appear explicitly in the GP equation (0) . These 
coupled equations (§)-(§) (along with expressions (Al) 
and (A2) in Appendix A for the collision integrals C12 
and C22) were derived in the semi-classical approxima- 
tion. They assume that the atoms in the thermal cloud 
are well-described by the single-particle Hartree-Fock 
spectrum ep(r,i) — p^ /2ni + U{r,t), where U{r,t) — 
C4xt(r) -I- 2g[nc{r,t) + n{r,t)]. However, they are ex- 
pected to contain the essential physics in trapped Bose- 
condensed gases at finite T, in both the collisionless and 
hydrodynamic domains. 

If the gas is weakly disturbed from equilibrium, one 
can consider linearized collective oscillations about equi- 
librium by writing 

$(r,t) = e-^'^-*/''[$o(r)+(5<i>(r,i)], (5) 

f{p,r,t) = f{p,r)+df{p,r,t). (6) 

Substituting (^ and (^) into Eqs. (|l)-(^, one can ob- 
tain the coupled dynamical equations for the damped os- 
cillations of the condensate 6^{r,t) and non-condensate 
(5/(p,r,t). We find that the damping term R{r,t) in 
(0) is finite even if we only keep the equilibrium part 
of /(p, r, i). This approximation simplifies the problem 
tremendously and should give a good first estimate since 
,5/(p,r,i)«/0(p,r). 

The thermal cloud influences the dynamics of the con- 
densate through two terms in (0) , the mean- field interac- 
tion potential 2gn{r,t) and the second-order colhsional 
term iR{r,t). Let us first consider the effect of the term 
2gh(r,t). In the static approximation, the term 2gno{r) 
gives rise to a temperature-dependent frequency shift of 
the condensate collective modes, but does not give rise 
to damping. Landau damping appears only when the 
dynamics of the non-condensate Sh(r, t) is accounted for 
(i.e. when we keep 5/(p, r, i)). An important point is 
that the collisional damping described by the second- 
order term ii?(r, t) in (pi) arises already within the static 
approximation (5/(p,r,t) = 0, in contrast to Landau 
damping. 



The damped collective mode equations for the con- 
densate can be obtained if we substitute the equilib- 
rium distribution /°(p,r) and density no(r) of the non- 
condensate into (g), giving us 



ih 



a$(r,i) 
dt 



2m 



- iRoir,t) $(r,i) 



V^ + C/cxt(r) + gndr, t) + 25^0 (r) 

(7) 



which describes the condensate motion coupled to the 
static thermal cloud. Here the damping term i?o(r, t) is 

fi 



Roir,t) 



dp 

2nc{r,t) J (2^3 



:Ci2[r(p,r),$(r,t)]. (8) 



Notice that Ro{r,t) now depends on time only through 
^{r,t), which has to be determined self-consistently by 
solving (|j). 

The equilibrium stationary solution of the coupled 
equations (||)-(||) is given by the solution of the gener- 
alized time-independent Gross-Pitaevskii equation 

iJo$o(r) = McO$o(r), (9) 

where we have defined the Hermitian operator 

Hn = -T^V^ + C/ext(r) + <?n,o(r) + 2giio{r). (10) 
2m 

Here, nco{r) — |<I>o(r)p is the static equilibrium density 
of the condensate and no(r) is the equilibrium density of 
the non-condensate. The equilibrium chemical potential 
of the condensate /ico is independent of position and can 
be written explicitly in terms of the equilibrium densities 
as 



A*cO — 



2m^nco(r) 



+ C^cxt(r) + gricoir) + 2gho{r), 



_ (11) 
where the first term on the right-hand side of (|I^) is the 
so-called quantum pressure of the static condensate wave 
function. 

To be consistent with the underlying kinetic model de- 
veloped in Ref. [|l|, the static thermal cloud is described 
using the single-particle HF spectrum (see above). 
Thus the equilibrium distribution /''(p,r) of the non- 
condensate atoms is given by the Bose-Einstein distribu- 
tion 



/°(p,r) = 



1 



3/3[pV2m+f7o(r)-Ao] _ 1 



(12) 



where (3 = l/ksT, fiQ is the equilibrium chemical po- 
tential of the non-condensate, and Uo{r) — f/cxt(r) -I- 
2g[nco{r) + no{r)]. As discussed in Ref. [Q, the static 
chemical potentials of the two components are equal 
(/io = Hco) in equilibrium. 

The equilibrium density of the non-condensate n.o(r) is 
obtained by integrating (n2h over momentum to give the 
usual result 



^oW = j3 93/2[Mr)]- 



(13) 



Here, A = (2nh'^P/my^^ is the thermal de Broghe wave- 
length and g3/2iz) is a Bose-Einstein function. The local 
fugacity is Zoir) = e'3('i°-'^o('-))_ 

We recall that in the TF approximation used in our 
earlier work ^ , (jl^) reduces to fico = Ucxt + gnco + ^gho 
and hence 20 (r) — e~^^"'''°^^\ When we keep the quan- 
tum pressure terms in (U^), the analysis is more compli- 
cated. The equilibrium solutions for the condensate and 
thermal cloud, given by (^ and (|l^), must be obtained 
self-consistently. For completeness, we outline our pro- 
cedure for a given value of the temperature T and with 
the total number of atoms N fixed (see also Rcf . ^^ ) : 

1. For a given value of the condensate population 
Nc < N, the equilibrium solution of the conden- 
sate is obtained by solving (^) numerically ||ll| . To 
start the procedure, we initially take Uq = 0. 

2. Using /Ltco and nco(r) — |$o(r)l'^ found from step 1, 
the non-condensate density no(r) is obtained from 
(|l§), using flo = ficO- If Mco > min{C/(r)}, which 
can occur for small condensates at temperatures 
close to Tbec, we set /lo = min{f7(r)}. 

3. Steps 1 and 2 are repeated until convergence is 
reached to the desired accuracy (in our calculations, 
when /ico stops varying up to an error of 10"^). 

4. The number of atoms in the thermal cloud is cal- 
culated N = J drno(r). We then repeat steps 1-3, 
varying Nc until the chosen total number of atoms, 
N = Nc + N , is obtained to the desired accuracy 
(in our calculations, we obtain N to an error of less 
than 10"''). 

The equilibrium values of nco(r), /^co = Aoi and no(r) 
obtained from the above procedure are used in the cal- 
culation of the damped collective oscillations of the con- 
densate discussed in the next section. 

Before proceeding, though, it is useful to emphasize 
several points about our approximate model, as described 
by (M) and (pi). Our emphasis in this paper, and the 
earlier TF version in Ref. Q, is on the calculation of 
damping due to collisional exchange of atoms between a 
dynamic condensate and a static thermal cloud. While 
we use the semi-classical HF excitations to calculate 
/°(p,r), which leads to ( [ll|) and (|lj), we do not ex- 
pect improved treatments of the static thermal cloud will 
greatly alter our estimates of the inter-component con- 
densate damping we are considering. There is a consid- 
erable literature on the calculation of the undamped con- 
densate mode frequencies at finite T. One finds that the 
frequency shifts are quite dependent on the specific ap- 
proximation y, Q, y used for the non-condensate mean 
field 2gnQ{r), as well as on other terms left out of (M) 
related to the anomalous correlation function rhoir) — 
('(/'(r)V-'(r)) We are not concerned with these questions 
here, but refer to Refs. Il2, 13 for further studies. 



III. BOGOLIUBOV EQUATIONS WITH 
DAMPING 



Using the explicit general expression for C12 given in 
( [A2| ) of Appendix A, it was shown in Ref. pi that Ro{r, t) 
can be simplified to (see also Ref. [[l4|) 



i?o(r,t) 



2ri2(r,i) 



[e-/3(Ao-ec(r.t))_;L], (14) 



where we have defined the relaxation rate due to inter- 
component C12 collisions 



1 



n2(r,t) (2 



7 /rfpi /rfp2 dpsdipc + Pi - P2 - Pa) 

X S{e, + ep,-ep,-ep,){l + f'l)f°f!^. (15) 



The condensate atom local energy is ec(r, t) = /ic(r, t) + 
^mv'^{r, t) with the non-equilibrium condensate chemical 
potential 



Aic(r,t) 



n'v^^^nAvT^ 



-|-[/cxt(r)+gnc(r, t)-|-2.gno(r) . 



2myjnc{v,t) 

(16) 
The condensate atom momentum is Pc — mwc, and 
fi — /"(i"> Pi)- We have introduced the usual condensate 
velocity de fined in t erms of the phase 9 of the condensate 
$(r,t) = v/nc(r,t)e''''(''^*) as v^ = hVe(r^t)/m. A closed 
set of equations for <i>(r, t) is given by ^ combined with 
(P) and (^ - ^. The damping term i?o(r, t) given in 
(14) vanishes when the condensate is in diffusive equilib- 
rium with the thermal cloud, when ec(r, t) = /i^o = fj-o- 

We obtain the linearized equation of motion for the 
condensate fluctuation 6^{r,t) by expanding all conden- 
sate variables appearing in (|^) to first order in S^{r,t). 
For simplicity we assume $o(r) is real, that is, Vco(r) = 0. 
Using (n), the condensate density nc{r,t) = |$o(r) + 
S^{r,t)f can be written nc(r, i) = nco(r) -I- Snc{r,t), 
where the linearized density fluctuation of the conden- 
sate is 



Sncir,t) = $o(i-)[(5$(r,t) +(5$*(r)]. 



(17) 



We can also simplify Ro{r,t) in (O) by noting that, to 
first order in the condensate fluctuations, we can write 
£c(r,t) — fJ.c{Tc,t) (neglecting the quadratic term (fac)^)- 
Using (O), the condensate chemical potential dlq) can 
be written as 



IJ.c{r,t) = HcQ + Snc{'r,t), 



(18) 



where the linearized fluctuation in the local condensate 
chemical potential is found to be given by 



Sfj.c{r, t) 



1 



[Hq — ^co) 



Snc{r, t) 



*o(r) 



g5n,{v,t). (19) 



where the operator Hi^ is defined in (pXJ). The first term 
in (M) arises from the dynamic quantum pressure of the 



condensate oscillation, while the second term comes from 
the mean field interaction. The latter is the only contri- 
bution kept in the TF approximation used in Ref. [g. 
Note that because the non-condensate is static, there is 
no contribution in (Km from the fluctuation in the HF 
mean field {2g Sh) of the thermal cloud. 

Utilizing the fact that ii^o — p-o in static equilibrium, 
Ro{r,t) in (|lj) can be simplified to 



6Roir,t) 



np 



2r?2(r) 



SHc{r,t), 



(20) 



where the "equilibrium" C12 collision rate is defined by 
(compare with (|l[ 



f 



2.9^ 



rUr) {2^fTJ 



dpi 



K 



v\ 



V2-V1 



2m 



dp2 / dp3S{pi - P2 - Ps) 
9nco){l + f?)f^fl (21) 



We can now obtain an equation of motion for the con- 
densate fluctuation 5$(r, t) by substituting (g) into (0) 
and using the above results 



ih 



dt 



£(5$(r,i)-|-gnco(r)(5$*(r,t) 

if [5$(r,t)+(5$*(r,t)], (22) 

^0 + 



We find it convenient to define the operators C 
guco - McO and 



2t°~ 



1 



(Hq - fie 



gnco 



(23) 



both of which arc Hermitian. The first term in (^3|) arises 
from the quantum pressure of the condensate fluctuation, 
ft is neglected in the Thomas- Fermi approximation, valid 
in the large N^ limit. The damping operator can also 
be written in an alternative form as F = {hf3 / Ati2)[C + 
gncoir)]- 

We now consider normal mode oscillations of the con- 
densate by expanding (5$(r, t) as 



5$(r,t) = Uj{r)e-''''/'' + v*{r)e'''^'/'', 



(24) 



where the excitation energy Ej and the collective mode 
amplitudes Uj{r) and Vj(r), in general, have real and 
imaginary components. Substituting this into ( |2^ ) and 
equating like powers of e^^"^^*/^, we obtain the damped 
finite-T Bogoliubov equations 

{C - ir)uj{r) + [gUco - iT)vj{r) = ejUj{r) 

{jC + ir)vj{r) + {gnco + ir)uj{r) ^ -ejVj{r). (25) 

The real and imaginary parts of Ej give the oscillation 
frequency loj = KeSj/U and the damping rate jj = 
—Imej/h, which vary with temperature T. These cou- 
pled equations are the main new result of this paper. 



They describe the damped collective modes of a conden- 
sate coupled to a static thermal cloud at finite T. ft is 
important to realize that the anti-Hermitian operator iT 
appearing in ( |25| ) has not simply been inserted "by hand" 
based on phcnomcnological arguments, as is done in the 
theories developed in Refs. ||l^, |l^, ^^, but has been 
derived explicitly starting from a microscopic (albeit ap- 
proximate) theory. The origin of this damping is the 
lack of collisional detailed balance between the conden- 
sate and non-condensate, which occurs when the system 
is disturbed from equilibrium. This mechanism is distinct 
from the Landau damping process, which arises from the 
mean field interaction between components when the dy- 
namics 5/ (p, r, t) of the non-condensate is taken into ac- 
count fill. 

In order to explore the properties of the coupled Bo- 
goliubov equations (25), it is useful to write it in matrix 
form 



Mwj 
where M = Mq + M' and 

Mo = 



= Eja^Wj, 

C gUco 
gUco jC 



M' 



-iT 



crs 



w'j(r) 



1 1 
-1 -1 



1 
-1 

Wj(r) 
^'i(r) 



(26) 



(27) 



(28) 



(29) 



(30) 



Since M' is not a Hermitian operator, the solutions of 
(Pq) are not in general orthogonal to one another and 
the energies Sj are expected to have an imaginary com- 
ponent describing damping. We can see this explicitly 
by deriving a generalized orthogonality relation for the 
solutions Wj ||l^, |9|. Multiplying (^6|) on the left-hand 
side by w] and integrating over position, we find 



dr'wj{r)Mwj{r) = Ej / drwl{r)a3Wj{r). (31) 



Taking the Hermitian conjugate of (|2^), replacing the 
index j by i, multiplying on the right-hand side by Wj, 
and integrating over position, we obtain 



drw'l{r)M^Wj{r) = e* / drw|(r)cr3Wj(r). (32) 



Subtracting (pll) from ( p2| ) gives us the desired orthogo- 
nality relation 



(e* -Ej) / drwl{r)a3Wj{r) = 2i / drwj{r)a3Twjir) 



(33) 



where we have used the Hermitian property Mq — M^ 
and used the fact that (M't - M') = 2zcr3f. 

If we neglect the damping term T ~ m (26) 



our description reduces to the usual coupled Bogoli- 
ubov equations at finite temperatures for the collective 
modes uqjM and voj (r) in the static Popov approxima- 

tion a 0, k M 



MqWqj = eajasWQj. 



(34) 



With r = in (|3^), we see that these solutions wqj obey 
the orthogonality condition 



is, 



Oi 



eoj) / drwl^a^woj = 0. 



(35) 



It then follows that the eigenvalues must be real {eoj — 
EqA, otherwise the solution woj of (|3J) would have zero 

norm ||2^. This property follows from the fact that Afg 
is a Hermitian operator. We note that at T = 0, where 
2ghQ = in M q, ('p4| ) reduces to the standard Bogoliubov 
equations |jl^, |19[. At finite temperatures, the simplest 
way of including the non-condensate is to use the temper- 
ature dependent condensate number N^iT) (but ignore 

the mean- field term 2577,0 (r) in Hq) H- As discussed in 
PI, the condensate normal mode frequencies at finite T 
are the same as at T = in the Thomas- Fermi limit 0, 
since the frequencies at T = do not depend on the value 
of Nc in this limit. Another effect comes from the mean 
field 2ghQ due to the non-condensate, which gives rise to 
a temperature-dependent energy shift to the excitation 
frequencies P, H|. However, as discussed at the end of 
Section II, a proper estimate of the shifts in the conden- 
sate normal mode frequencies requires a better model. 

According to (33), the inclusion of the anti- Hermitian 
damping operator iV in the normal mode equations in- 
troduces an explicit imaginary part in the eigenencrgies 
Sj, without requiring the solutions to have zero norm. 
However, these damped solutions Wi are not required to 
be orthogonal, 
given by 



From (p3D, the imaginary part of Ej is 



Im£i 



- J drwjasTwj 
J drWjffsWj 



(36) 



If the operator F in (g3|) was a scalar position- 
independent constant F = 7, we can take it out of the 
integral in (p3), and in this case the orthogonality con- 
dition for Wj is recovered. 



(s:-e,-2^,)/d™t.3^,.0. 



(37) 



This states that for i — j, the imaginary part of each of 
the excitation energies is Ini£j = —7 (unless the norm of 
Wj is zero). For i ^ j, ( pTI ) shows that the eigenmodes 
Wi and Wj are orthogonal as long as Ree^ ^ ReSj. 

Since the damping effect of C12 collisions on the con- 
densate oscillations is weak, we can treat the damping 



term M' in ( |28| ) using first order perturbation theory. 
More precisely, we assume 



fiujj > 7j , 



(38) 



where 7-,- = — Irne,- is the damping rate due to C12 colli- 
sions given by ( |36| ) . Proceeding as usual in perturbation 
theory, we introduce a small expansion parameter 



M = Mq + AM', 
and expand the energies and eigenmodes as 



Wi 



Eoj + Aeij + A"e2j + ■ • 

Wqj + Xwij + \^W2j + 



(39) 



(40) 
(41) 



Substituting ( p9|) - (^) into (|2g) and equating like powers 
of A, we obtain a hierarchy of equations. Note that since 
M' is also radially symmetric, it does not break the [21 + 
l)-fold degeneracy of the magnetic sublevels. The zeroth 
order equation simply reproduces ( |34| ) for the undamped 
modes wqj • The first order correction is obtained from 



(A/o - eoja:i)wij = (eijcra - A/') 



Wo J. 



(42) 



If we multiply (^) on the left by wJ ■ and integrate over 
position, the left-hand side of ( p2[ ) vanishes using (p3). 
This leaves us with the following result for the perturbed 
energy shift. 



Eij 



J drwLM' 



6j^^^ ^oj 



j drwl-a^WQj 



(43) 



Taking the Hermitian conjugate of (^), multiplying on 
the right by woj and integrating over position, we obtain 
the analogous expression for the complex conjugate eJ 



J drwhM''^WQj 



(44) 



From 

first order) 



Jdrwl^aswoj 
and (PI), we can obtain tlie damping (to 



loj = -Imeij = 



J drwljastwQj 
J drwl^aswoj 



(45) 



Here we have used (M'^ - M') = 2io-3f. The result in 
( [45| ) is consistent with expression ( p6|) , had we simply ex- 
panded ^a) about wqj. We can also obtain the real part 



of £ij from (|4^) and ([44| ) and we find that this vanishes 



(see Appendix B), i.e. 



Reeij = 0. 



In summary, first order perturbation theory gives 



eoj -ijor 



(46) 



(47) 



Using these results, the expansion given in {UM for the 
condensate oscillations can be written in the more ex- 
plicit form 



(5$(r,i) 



— g-70jt/?i| 



uoj{r)t 



-ieojt/h 



^Oj 



(r)e> 



eojt/h) 



(48) 

Using expression (|45| ) to calculate 7oj, one can check that 
the criterion in (|38D is well satisfied. The difference be- 
tween the pcrturbative results and the exact numerical 
solution of ( p5| ) is always less than 1 % in our calculations. 
The results of the present calculation can be compared 
to those in Ref. |^], where we used a simplified model to 
calculate the damping rates jj. We used the Thomas- 
Fermi approximation (neglecting the kinetic energy pres- 
sure in the solution of $0 and Hco) and we also neglected 
the contribution 2gnQ to the mean field interaction. We 
also neglected the kinetic energy pressure in the solution 
of the normal modes, which allowed us to use the T — 
Stringari normal mode solutions H as a basis with which 
to treat the effect of SRq to first order. The main result 
of Ref. H was the expression for the damping rate 



7j 



hJdrSnlir)/T'ir) 
2 JdrSnUr) ' 



(49) 



Here Snsj{r) is the density ffuctuation associated with 
a Stringari normal mode [^ |^ at finite T and 1/r' 



{gnM/kBT){l/T^2)' with I/t^^ given in ^. This damp- 
ing rate is easy to evaluate since the equilibrium form for 
the condensate $o(r), as well as the Stringari normal 
modes Sngj (r) , are both given by simple analytic func- 
tions. 

The TF result (^9|) can, of course, be obtained directly 
from (Ea) by taking the Thomas-Fermi limit. In this 
limit, the first term in (E3^ for T. which describes the 



quantum pressure of the oscillation, can be neglected. 
We are left with r(r) = ?i/2r'(r), and hence ( [45| ) reduces 
to 



7oj 



h Jdrwlja3Woj/T'{v) 
2 J drwLasWQ 



(50) 



In Ref. ||2l|, the product wl ■(r)cr3Woj(i") — |Mj(r)|^ 
|wj(r)p is shown to satisfy the following relationship 



K(r)p-K-(r)p=Re: 



-2im 



5(j)*{v)5nj 



(51) 



where 5nj{r) is the density fluctuation of mode j, and 
d(j)j{r) is the corresponding fluctuation in the velocity 
potential. In the Thomas-Fermi limit, the finite T fluc- 
tuations in the density and velocity potential still have 
the simple relationship |23] 



*s,(r) = 



-w 



muJsj 



Snsjir), 



(52) 



where wgj is the Stringari frequency of mode j. Putting 
these results together, we have 



wlj{r)a3Woj{r) = ^Snlj{r) 



(53) 



and hence (|50| ) reduces to ([1 91). 

In the next section, we compare the damping given 
by ( p^ ) to the result given by (^^ for a trapped gas. 
It is useful to apply ( p5| ) to the case of a homogeneous 
gas, where the modes are Uk{r) — Ufec''''" and 'yk(i") = 
tifee*'''" and the energies are given by the usual Bogoliubov 
expression 
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(54) 



Usin g th e expression for F given in 
rate (Eg) reduces to 



231), the damping 
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tip 
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Am 



guco , (55) 



where T12 , as given in (p^) , is independent of position for 
a uniform Bose gas. The first term proportional to k^ 
comes from the quantum pressure (i.e. the first term in 
(p3|)), which is omitted in the Thomas- Fermi result given 
in Ref. §. 
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FIG. 1: Condensate fraction versus temperature for A'^ — 
2 X lO" atoms. The solid line corresponds to the method 
described in steps 1-4 in Section II. The Thomas- Fermi ap- 
proximation is used to calculate the dashed line, as described 
in Ref . H, bs] . The dot-dashed line is for an ideal Bose gas in 
the thermodynamic limit. 
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FIG. 2; Condensate mode frequencies versus temperature for 
N — 2 X 10® atoms. The dashed lines are the T = modes 
but for a changing number of condensate atoms Nc{T). The 
solid line is the real part of the eigenenergies £i obtained by 
solving (P6|). 



IV. NUMERICAL RESULTS 

We now turn to an explicit calculation of the inter- 
component damping of collective modes for a dilute Bose 
gas in a spherical harmonic trap f7ext(r) — ^rnuj^r'^ . We 
choose the same set of physical parameters as given in 
Ref. [||: the frequency of the trap is ujo/2tt — 10 Hz, the 
scattering length for ^^Rb is a = 5.7 nm, and the total 
number of atoms is iV = 2 x 10^. 



A. Equilibrium solution 

In Fig. 1 we plot the condensate fraction versus tem- 
perature for the case of A^ = 2 x 10^ atoms. The solid line 
corresponds to the full calculation of Nc{T) as described 
in steps 1-4 of Section II. The dashed line is obtained 
using the Thomas-Fermi approximation for $0 and ne- 
glecting the effect of 2gno in Hq and U{r), as described 
in Refs. H, E3|. The two curves deviate only slightly as T 
approaches Tbec- This is consistent with previous stud- 
ies [^ |24 1 , which found that the thermodynamic quanti- 
ties of a dilute Bose gas are quite insensitive to the level 
of approximation used to treat the non-condensate spec- 
trum. In Fig. 1, for comparison, we also show the result 
for the case of the ideal gas in the thermodynamic limit, 
Nc/N = [1 - {kBT/hcjo)^] (t), given by the dot-dashed 
line. When interactions are included, the effective con- 
densation temperature Tbec is slightly lower Q. 



B. Collective mode frequencies 

We plot the collective excitation frequencies uji in 
Fig. 2 for the breathing {n = 1,Z = 0, m = 0}, dipole 



{0,1, to}, and quadrupole {0,2,to} modes. The solid 
line corresponds to the numerical solution of the full cou- 
pled Bogoliubov equations in (Eq). For comparison, the 
dashed lines are obtained from the solution of ( p4[) set- 
ting 2gno = and using the value of Nc{T) (as given by 
the dashed line in Fig. 1). In the Thomas- Fermi Stringari 
limit, the frequencies are given by the T = results, at 
all temperatures. 

As T approaches Tbec, Nc{T) becomes very small 
and the solutions given by the dashed lines go over to 
the harmonic oscillator eigenstates — the breathing and 
quadrupole modes become degenerate at uJi/cJn = 2. As 
T decreases, the number of atoms in the condensate in- 
creases, so that the dashed line approaches the T = 
Stringari frequencies. The dipole mode frequency differs 
from the value Ui = luq expected for the Kohn mode since 
the thermal cloud is treated statically (see also Ref. [Q). 

The deviation of the full solution (given by the solid 
lines) from the dashed lines is due entirely to the mean 
field of the static thermal cloud 2gfiQ. The upward 
shift in the frequencies (as illustrated in Fig. 2 for 
iV = 2 X 10^) increases with A''. Early studies [|,|, |] 
of collective modes in the static Popov approximation 
considered much smaller systems, N ~ lO'^. In more 
recent work |I3], the mode frequencies were computed 
for a much larger system of iV = 2 x 10^ atoms. Al- 
though the thermodynamic quantities are very insensi- 
tive to the level of approximation used to treat the spec- 
trum of the non-condensate atoms Jl^, |2J], the lowest- 
lying collective mode frequencies arc clearly more sensi- 
tive. While we feel it is useful to show what our simple 
model gives for the frequency shifts, the shifts found in 
Refs.|g] and [||, based on calculating no(r) more sclf- 
consistently, are much smaller than our model predicts. 
A good estimate clearly requires a more realistic the- 
ory O, O, as discussed at the end of Section II. 
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FIG. 3: Damping rates versus temperature for N — 2 x lO'' 
atoms. In the main plot, the solid line corresponds to the full 
calculation of "/i — —Irasi for the breathing mode as obtained 
by solving (Efs). The simplified TF expression (Ii9|) is given 
by the dashed line. In the inset, the full calculation of ji 
is shown for the breathing (solid), dipole (dot-dashed), and 
quadrupole (dashed) modes. 



C. Damping due to inter-component collisions 



In Fig. 3, we plot the damping rates ji for N = 2x 10^ 
atoms. In the main graph the damping of the breath- 
ing mode is shown, where the solid line comes from a 
full calculation based on (|2^) and the dashed line is cal- 
culated using the simplified expression in ([49| ) derived in 
the Thomas-Fermi limit H . As a reference point, we have 
indicated the temperature at which the condensate pop- 
ulation Nc(T) reaches 10^ atoms, where the ratio in (l[) is 
rno/RiF = 0.33. The largest deviation, as expected, oc- 
curs close to Tbec where the Thomas-Fermi approxima- 
tion (|^) used in deriving (^ ) starts to break down. The 
damping becomes very large as Nc goes to zero close to 
the condensation temperature. This increase arises from 
the kinetic quantum pressure of the collective oscillation, 
i.e. the first term in (|2S 



In the inset of Fig. 3 we graph the damping rates for 
all three modes considered in Fig. 2. The similarity em- 
phasizes the fact that the damping is fairly insensitive to 
the detailed form of the condensate normal modes, which 
suggests that the inter-component damping shown in Fig. 
3 will not be significantly modified by a more accurate 
description of the static thermal cloud, in contrast to the 
frequency shifts. 

Clearly as we approach TbeCj the condensate is becom- 
ing small while the thermal cloud is becoming more dom- 
inant. In this transition region, the nature of the excita- 
tions in a trapped Bose gas is quite complex. Since the 
condensate is disappearing, so are the collective modes 
and a theory based on (0) becomes inadequate. 



V. DISCUSSION AND OUTLOOK 

In summary, we have derived damped collective mode 
equations (|2^) for a condensate interacting with a static 
thermal cloud starting from the generalized GP equa- 
tion derived in Rcf. m. This mechanism is due to the 
lack of collisional detailed balance between the two com- 
ponents when the condensate is disturbed from equilib- 
rium. In a previous article [p|, we evaluated this damp- 
ing in the large Nc Thomas-Fermi limit and neglected 
the mean field of the non-condensate 2gho. In Section 
IV we presented results of an explicit calculation of ( [25| ) 
and verified that the simplified TF expression ( p^ ) agrees 
quantitatively with the damping obtained from the full 
calculation, except very near Tbec where Nc{T) is be- 
coming small. 

Starting from the finite T generalized GP equation (1^), 
we formulated our calculation in terms of coupled Bo- 
goliubov equations (]25|), in which the inter-component 
collision damping arises explicitly. We have used this 
opportunity to discuss some of the formal properties of 
these equations when damping is present. Such equations 
may be of interest in other problems when dealing with 
damping from the thermal cloud. 

The effect of damping in the time-dependent Gross- 
Pitaevskii equation has been discussed in several differ- 
ent contexts |l5| |l^, |l^. Many years ago, Pitaevskii |1^ 
developed a phenomcnological model for superfluid he- 
lium to describe the evolution of the superfluid toward 
equilibrium in the two-fluid hydrodynamic regime. This 
model has points of contact with the work in Rcf. M . In 
Ref . [0| , a similar scheme was used to discuss damping of 
condensate modes in the collisionless region. In Ref. p6J , 
damping terms were introduced into the GP equation to 
account for output coupling of condensate atoms from 
the trap as well as the exchange of atoms between the 
condensate and thermal cloud. Finally, in quite a differ- 
ent context, damped collective modes were calculated in 
Ref. [g5| for a condensate in an optical cavity. 

Taking into account the dynamic mean fleld coupling 
between components, the collective modes of the system 
become a hybridization of condensate oscillations and 
non-condensate oscillations, so that one gets essentially 
a pair of in-phase and out-of-phase oscillations for each 
collective mode. For a given mode symmetry, the out- 
of-phase mode consists mostly of condensate oscillations 
with the collective motion of the non-condensate being 
less significant [|l[. Our model calculations (and those 
of Ref. H ) should be quite good for the intercomponcnt 
Ci2 damping of such out-of-phase modes. In the other 
extreme, the in-phase Kohn mode pi involves both the 
condensate and non-condensate moving together in local 
equilibrium, in which case -R(r, t) = 0. 

The dynamic coupling between the condensate and 
thermal cloud gives rise to damping of the condensate 
oscillations (known as Beliaev and Landau damping), 
which are quite distinct mechanisms from that consid- 
ered in the present paper. Beliaev and Landau damping 



arise from dynamic mean- field effects, rather than colh- 
sion terms such as C12 or C22 (for further discussion and 
references, see Refs. g |lj). In a complete theory the 
condensate collective mode damping would be given by 
7tot = 7l + 7s + 7Cj that is, Landau, Beliaev, and col- 
lisional (due to C12) damping, respectively. Our model 
gives a good estimate of the C12 damping 7c. In our ear- 
lier work [g| where we used the TFA, we compared the 
C12 damping directly with Landau damping and found 
that 7c /7l ~ 0-3 to 0.5 for typical trap parameters. 

We thank Sandy Fetter, Tetsuro Nikuni, Reinhold 
Walser, and Milena Imamovic-Tomasovic for useful dis- 
cussions. This work was supported by NSERC. 



From (Wa) and (44) we can write 



Ree 



ij 



2 Jdrwl^aaWoj 



0-2 = 



-1 

1 



(Bl) 



From the explicit form of M' given in (Jlg), we obtain 
(M' -^ M't) = 2ia2t, where 



(B2) 



APPENDIX A: COLLISION INTEGRALS 

The two collision terms on the right-hand side of 
are given by M 



C22[/] = -^ Jdp2 Jdps Jdp, 



(2 
X S{p + P2 - P3 - P4)^(ep + £p2 - £p3 - ^Pi) 

X [(1 + /)(! + /2)/3/4 - ,//2(l + /3)(1 + fi)] , 

(Al) 

which is the usual Uehling-Uhlenbeck form of the Boltz- 
mann collision integral describing collisions between ther- 
mal atoms, and 



Cl2[f,^ = 



2g'^nc 



dpi dp2 dp. 



X S{mVc + Pl - P2 - P3)^(ec + £pi - £p2 - Epa) 

X [<S(P-Pl) -(^(P-P2)-^(P-P3)] 

X [(l + /l)/2/3-/l(l + /2)(l + /3)], (A2) 

describing the coUisional exchange of particles between 
the condensate and non-condensate. 



APPENDIX B: PROOF OF EQ. 46 



Then (Bl) becomes 



Ree 



ij 



i J drwlja2TwQj 
j dvwl^a^WQj 



(B3) 



We next define w' — Vwqj , which we expand in the basis 
of the undamped eigenmodes 



"'-E^' 



WQi 



(B4) 



Substituting this into ( B3 ) gives 



Reeij — i y^ ' 



J drwlja2Woi 
J drwl^a-iWoj 



(B5) 



We make use of the general result |I8| 



drwlj{r)a2Woi{r) ^ j dr[v*{r)u.,{r) ~ u*{r)vt{r)] 

(B6) 



We give a formal proof that there is no first-order cor- 
rection to the condensate frequencies from F in (25) 



which holds for all i and j. Using this in (B5) gives the 
final resuh (Mj. 
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